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Abstract. In recent years there have been new insights into the integrability of 
quadrilateral lattice equations, i.e. partial difference equations which are the natural 
discrete analogues of integrable partial differential equations in 1+1 dimensions. In the 
scalar (i.e. single-field) case there now exist classification results by Adler, Bobenko 
and Suris (ABS) leading to some new examples in addition to the lattice equations "of 
KdV type" that were known since the late 1970s and early 1980s. In this paper we 
review the construction of soliton solutions for the KdV type lattice equations and use 
those results to construct A^-soliton solutions for all lattice equations in the ABS list 
except for the elliptic case of Q4, which is left to a separate treatment. 



1. Introduction 

The study of integrable partial difference equations (P AEs) dates back to the pioneering 
work of Ablowitz and Ladik, [1] , and of Hirota, [18] , motivated partly by the search for 
numerical finite-difference schemes that are discrete in time as well as in space. Such 
systems constitute discrete analogues of soliton type PDEs. In subsequent work the 
Lie-algebraic approach of the Kyoto school, [12], on the one hand, and the approach 
using a so-called direct linearisation method, cf. [221 EZ], on the other hand led to new 
systematic constructions for such systems. In recent years they have been reinvestigated 
from various points of view, including reductions to integrable dynamical mappings 
[251 [26], cind associated finite-gap solutions [10, [lH [IS] , reductions to discrete Painleve 
equations, [211 [23] and similarity reductions, whilst soliton solutions arose as a direct 
corollary from the original constructions mentioned above. 

The property of multidimensional consistency, first set out explicitly in pi], cf. 
also [11], lies implicitly within lattice equations which have the interpretation of a 
superposition principle for Backlund transformations. This was identified in [23] as 
the property constituting the precise discrete analogue of the existence of hierarchies of 
nonlinear evolution equations, and hence of integrability. The property has subsequently 
been used by Adler, Bobenko and Suris [3l [1] as a classifying property. Within certain 
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additional conditions they produced a full list of scalar quadrilateral lattice equations 
which are multidimensionally consistent. This list, which surprisingly is quite short, is 
reminiscent of Painleve's list of transcendental equations in the case of second order 
ODEs possessing the property of non-movable singularities of the general solution. We 
reproduce this list from [3] belo\Jil containing three groups of equations: the Q-list the 
H-list and the A-list: 



Q-list: 



Ql : p(u — u){u — u) — q{u — u){u — u) = S'^pq{q — p) (1.1a) 
Q2 : p{u — u){u — u) — q{u — u)(u — u)+ pq{p — q){u + u + u + u) = 

= pq{p - q){'p - pq + q^) (1.1b) 
Q3 : p(l — g^)(nn + uu) — g(l — p^){uu + uu) = 

= if - q') [{nu + n§) + ^) (1-lc) 



Q4 : p{uu + uu) — q{uu + uu 

o o 

pQ - qP 



1 °9 °9 

1 — p^g^ 

where = p^ — + 1 , = — + 1 . 



(^(uu + uu) — pq{l + uuuu)^ (l-ld) 



H-list: 



HI: {u — u){u — u) = p — q (1.2a) 
H2 : {u-u)(u-u) = {p-q){u + u + u + u) +f -q^ (1.2b) 
H3 : p{uu + uu) — q{uu + uu) = 6{q^ — p^) (1.2c) 



A-list: 

Al : ]9(n + n)(n + ^I) — g(ti + ?I)(?i + 2^) = 5^]5g(p^ — g^) (1.3a) 
A2 : p{l — q'^){uu + uu) — q{l — p'^){uu + uu) + {p'^ — q'^) (1.3b) 



The notation we have adopted here and in earlier papers is the following: the vertices 
along an elementary plaquette on a rectangular lattice contain the dependent variables: 

U := "Un.m+l; U = '^n+l.m+l) 

which schematically are indicated in Figure [H in which p and q denote lattice parameters 
associated with the directions in the lattice (measuring the grid size in these directions): 

I Note that the Q4 equation as given by [3] , which was first found by Adler in [2] , is different from the 
one in (jl.ldp , which in tliis form was first presented in |16j . 
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Figure 1: Arrangement of the shifted dependent variable on the vertices of a 
quadrilateral and association of the lattice parameters to the edges. 



o o 

u ^ u, u ^ u. The lattice parameters associated to the two directions on the 
lattice play a central role in the notion of multidimensional consistency. Specifically 
they parametrise the family of equations which are compatible on the multidimensional 
lattice. The equations in the Q-list are related by degeneration (of the elliptic curve 
associated to Q4) through the coalescence scheme illustrated in figure [2J The H-list and 



Q35=i 




Q35=o 





Q4 




Q2 




Ql5=l 




Ql5=0 







Figure 2: Coalescence diagram for equations in the ABS Q list (this was first given 
by Adler and Suris in [5]). 

A-list appear also by degeneration from the Q-list, leading in principle to an extension 
of Figure [21 and we will make some of these connections explicit in the later part of the 
present paper. Connections between the ABS equations are not limited to coalescence 
by degeneration, there also exist Miura and Backlund type transformations connecting 
distinct equations in the list, some transformations of Backlund type were discussed 
recently in [9]. 

The present paper is part I of a sequence of papers dedicated to closed-form A^- 
soliton solutions of the lattice equations mentioned above. Part II, by Hietarinta and 
Zhang, [TTj, will be dedicated to the Casorati form of the same solutions, establishing a 
different methodology. The results in the present paper are based on a Cauchy matrix 
structure which goes back to earlier work [221 EZj, on the so-called direct linearization 
approach. In section 2 we review the application of that approach to soliton solutions of 
lattice equations "of KdV type" , which were cases already known in the earlier papers 
mentioned. Specifically these include the lattice potential KdV, potential modified KdV, 
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Schwarzian KdV and an interpolating equation between them which has been referred 
to as the NQC equation (cf [2H])- (These equations are equivalent to the equations 
HI, H35=o, Ql<5=o and Q35=o respectively up to gauge transformations.) The machinery 
introduced in section 2 will comprise several objects, and relations between them, which 
we will need throughout the remainder of the paper, because the iV-soliton solution for 
the full Q3 equation can be expressed in terms of those quantities. A constructive proof 
of this A^-soliton solution is presented in section 3, whilst in section 4 we show that 
this A^-soliton solution is consistent with the Backlund transformation (i.e., defining the 
latter by a copy of the lattice equation itself we establish the precise relation between the 
A^- and A^+l-soliton solution in terms of this BT). As an upshot of the present treatment, 
it becomes apparent that the natural parametrisation of Q3 turns out to involve already 
an elliptic curve, whose branch points can be viewed as lattice parameters associated 
with additional lattice directions. Degeneration from Q3 in this parametrisation by 
bringing together one or more branch points of the curve yields one or other of the 
equations on the ABS list. In section 5 the degenerations from this parametrisation of 
Q3 are given in detail and we construct the A^-soliton solutions for all the equations 
in the ABS list, except for the elliptic case of Q4. This reveals, on the level of the 
A^-soliton solutions, a deep connection between the equations expressed in terms of the 
basic quantities introduced in section 2. Section 6 contains a discussion of the results 
and concluding remarks. 



2. A^-Soliton solutions of KdV Type Lattice equations 

Among the equations in the ABS list we distinguish a subclass which we call of KdV 
type, and they comprise the lattice equations which have appeared already many years 
ago in the context of discretzations of the KdV equations and its counterparts, cf e.g. 
fL8 \ [22l [27] . In this section we introduce objects from which we may construct solutions 
of the KdV type lattice equations. 



2.1. Cauchy matrix and recurrence structure 

We will start by introducing the following Cauchy type matrix: 

which will form the core of the structure which we will develop. In (12. ip the q, kt, 
[i = 1, . . . , A^) denote two sets of A^ nonvanishing parameters, which we may chose 
freely (apart that we have to assume that ki + kj 7^ 0, Vz, j = 1, . . . , A^, in order to avoid 
difficulties with the numerators in the matrix M). These parameters are assumed not 
to depend on the lattice variables, i.e. on the discrete independent variables n and m. 
The dependence on the latter are incorporated wholly in the functions pi, the so-called 
plane wave factors, which have the form: 
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where the p° are initial values, i.e. constant wih regard to the variables n, 

Let us now introduce for convenience the following notation. Let K denote the 
diagonal N x N matrix containing the parameters ki on the diagonal, and introduce a 
column vector r, containing the entries pi, and a row vector *c, containing the entries 
Ci, i.e. 



K 



( k. 



\ 



kN / 



/ Pi \ 

P2 



r 



Cat J 



:2.3) 



\PN / 

It is easily checked that we have from the definition (12. ip immediately the relation: 

MK + KM = r^c, (2.4) 



where significantly the dyadic on the right hand side is a matrix of rank 1. 

Next we establish the dynamics in terms of the matrix A4", which follows from the 
definition (12.11) together with the dynamical equations for pi (12.21) which are simply 



~ m ^ P + ki 

Pi = TpPi = pi[n + 1, m) = — Pi 

p — Ki 



Pi 



Tqpi = pi{n,m + 1) 



q + kj 
q-ki 



P., (2.5) 



where Tp, Tg denote the elementary shift operators in the lattice in the directions 
associated with lattice parameters p and q respectively, T_p = and T_g = 
denoting their inverses. A straightforward calculation then yields the following relations 



M {pi + K) - {pi + K) M = r , (2.6a) 
M {ql + K)-{ql + K)M = 9^c, (2.6b) 

in which 1 is the N x N unit matrix, and where we used the obvious notation that the 
shifts ~ , ^ act on all the relevant objects depending on n, m by the respective shifts 
by one unit in these independent variables. In addition, we have the adjoint relations 

{pi - K)M - M {pi - K) = r^c , (2.6c) 
{ql - K)M - M {ql - K) = r . (2.6d) 

§ Importantly, the can still contain an in principle arbitrary number of additional discrete 
exponential factors of the form given in (|2.2p . i.e. depending on additional lattice variables, each 
of which would be associated with its own lattice parameter. Thus, a more general form for the pi 
would be 

n (^)^ 

containing an arbitrary number of lattice variables Ui, each with lattice parameter pi, labelled by some 
index ly. Every statement derived below can be extended without restriction to involve any choice 
of these variables, and in particular this implies the multidimensional consistency of the equations 
derived from the scheme. Thus, wc will say that these results can be covariantly extended to the 
multidimensional case. 
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The equations (I2.6p encode all the information on the dynamics of the matrix M, w.r.t. 
the discrete variables n, m, in addition to (12 ■4p which can be thought of as the defining 
property of M. 

Now we introduce several objects involving the matrix M, in terms of which we 
can define the basic variables which will solve the relevant lattice equations. Thus, we 
introduce the determinant 

/ = /n,n^ = det(l + M) , (2.7) 

which we will identify later as the relevant r-function (obeying Hirota-type bilinear 
equations), as well as the following quantified 

w» = + K'r (2.8a) 

= 'cK^ (1 + My^ (2.8b) 
Siid) = t^x^ (1 _/v^)-i r , (2.8c) 

for 2, j G Z (assuming that none of the parameters fcj is zero). Thus, we obtain an 
infinite sequence of column vectors n*^*\ of row vectors 'w^-'-' and a infinte by infinite 
array of scalar quantities S^'^'^\ An important property of the latter objects, which can 
also be written as 

^(^j) = t^K' ^{i) = t^U) X' r , (2.9) 
is that they are symmetric w.r.t. the interchange of the indices, i.e. 

Sii,3) = s^^''\ (2.10) 

provided that the constants Ci, p'^, k = 1 . . . N, are all nonzero. 

We shall now derive, starting from (12. 6p a system of recurrence relations which 
describe the dynamics for the quantities defined in (12.81) . Once this recursive structure 
is established, we will single out specific combinations of the S^'^'^'' in terms of which we 
can derive closed form discrete equations. In fact, from the definition (]2.8ap . using the 
fact that 1^ is a diagonal matrix, we have: 

K' r = {1 + M) u^'^ K'r = {l + M)u^^ 

pi- K ^ ' 
K' {pl + K)r = {pi- K){l + M)u'''^ 
= {pi - K) u^^ + {pi -K)M u^^ 
= {pi - K) u^'^ + M{pl- K) u^'^ + r *cu^^ 



II Here and earlier we use the symbols *c and *m to denote adjoint row vectors, which in the case of 
the first is just the transposed of the column vector c, i.e. *c = c^, but in the case of *m is not simply 
the transposed of the column vector u, but rather a quantity that is defined in (|2.8bp in its own right. 
Thus, the left super-index "t" should no be confused with the operation of transposition, but rather 
indicate a new object obeying some linear equations associated with the vector u. 
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where in the last step use has been made of fl2.6cp . Using now fl2.8cp we we conclude 
that 

pK^ r + K^+^ r = {l + M){pl-K) u^''^ + S^^'^V . 

Multiplying both sides by the inverse matrix (1 + M)~^ and identifying the terms on 
the left hand side using (f l2.8ap . we thus obtain 

(pi - K) u^'^ = (1 + My^ \pK' r + r - ^(^'°V 

= + - ^(^-o^n^o) . (2.11) 

Thus, we have obtained a linear recursion relation between the objects w'-*^ with the 
objects S'*-*'-'^ acting as coefficients. In quite a similar fashion we can derive the relation 

{pi + K) n« = pu^"^ - u^'+^^ + 5(^'0)n(°) , (2.12) 



in fact by making use of (12.6bp in this case. Eq. fl2.12p can be thought of as an inverse 



relation to (12. lip , noting that the ~-shifted objects are now at the right hand side of 
the equation. Multiplying both sides of either (12. lip or (12.120 from the left by the 
row vector '^cK^ and identifying the resulting terms through (I2.8cp . we obtain now a 
relation purely in terms of the objects S^'^'^\ namely: 

'cK' {pi - K) u^"^ = 'cK' \pu^'^ + u^'+^^ - 

^ pS^''^^ - = pS'^''^'> + - (2.13) 

using the fact that S*^*'"-* is just a scalar factor which can be moved to the left of the matrix 
multiplication. Thus, we have now obtained a nonlinear recursion relation between the 
S^"^'^^ and its ~-shifted counterparts. 

In a similar fashion, multiplying eq. (I2.12p by the row vector *c we obtain the 
complementary relation: 

pSiij) + = psi^j) _ + 5(i,o)^(o,i) ^ ^2.14) 

however this relation can be obtained from the previous relation (I2.13P by interchanging 
the indices, using the symmetry (I2.10p . 

Remark: Combining both eqs. (12. lip and (12.120 . using also (I2.14p . the following 
algebraic (i.e. not involving lattice shifts) recurrence relation can be derived: 

K^wW = 1.(^+2) ^ ^(^,1) ^(0) _ ^(^,0) ^(1) ^ (2.15) 
which in turn gives rise to the following algebraic recurrence for the objects S^^'^\ namely 

5'(i,i+2) ^ _^ g{i,l) g{0,j) _ g{i,0) _ (2.16) 
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To summarise the structure obtained, we note that starting from the Cauchy matrix 
M defined in (12. ip . depending dynamically on the lattice variables through the plane- 
wave factors Pi given in f l2.2p . we have defined an infinite set of objects, namely column 
and row vectors n*^*^ and *u^''\ and a doubly infinite sequence of scalar functions 5"^''-''', 
all related through a system of dynamical (since it involves lattice shift) recurrence 
relations. Obviously, all relations that we have derived for the ~-shifts (involving lattice 
parameter p and lattice variable n) hold also for the ^-shifts, simply by replacing p by 
q and interchanging the roles of n and m. Thus, for the scalar objects S^^'^^ we have the 
following set of coupled recurrence relations: 

pSi^j) _ s(',3+i) = psi^,j) + _ ^{i,o)^{o,i) ^ ^2.17a) 

pSi^'i) + = _ + 5(i.o)^(oj) ^ {2A7h) 

qS(''^) _ si',j+i) = ^^(^j) + _ ^{i,o)^{Oj) ^ ^2.17c) 

+ ^(ij+i) = qs'^^,J) _ + 5{i,o)^{Oj) _ {2A7d) 

We proceed now by deriving closed-form lattice equations for individual elements chosen 
from the S^'^'^^ as functions of the variables n, m. 

2.2. Closed Form Lattice Equations 

We start with the variable S^^'^\ for which we can derive a partial difference equation 
as follows. In fact, subtracting fl2.17cp from fl2.17al) we obtain 

pS^''^^ - qS^'''^ - S^''^+^^ + S^'''+^^ = {p- q)S^'''^ - (5(*'°) - . (2.18) 



On the other hand, taking the ^-shift of fl2.17bp . let us refer to it as fl2.17bp . and 



subtracting from it the -shift of fl2.17dp . i.e. f l2.17dp . we obtain 

_ _ ^ _ ^ ^ ^ ^{O.j) 

pS^'^^^ - qS^''^^ + ^(*'^'+^) - ^(^■^■+^) = (p - g)5(,j) + (^(^'°) - S^''^^)S . (2.19) 

Combining both equations, the terms which have a shift in their second index drop out 
and we obtain the equation: 

{p + q){S^''^^ - S^''^^) = {p- q){S^''^^ - f'''\ + (^('''"^ - S^(''°))(5(°'^') - . (2.20) 

Setting now i = j = in the last formula, we see that we get a closed form equations 
in terms of w = S^^''^\ This yields, after some trivial algebra the equation: 

{p + q + w — w){p — q + w — w) = p"^ — q'^ , (2-21) 

which is the lattice potential KdV equation, which has appeared in the literature in 
various guises, cf . [HI [22], [27] , notably as the permutability condition of the Backlund 
transformations for the KdV equation, cf. [29]. Curiously, this integrable partial 
difference equation can also be traced back to numerical analysis, where it has appeared 
in the form of the e-algorithm of Wynn, [31], as an efficient convergence accelerator 
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algorithm. In the present context of the structures arising from the Cauchy matrix we 
have estabhshed here an infinite family of solutions of the form 



w 



= *c(l + M)"ir (2.22) 



which constitute the A^-soliton solutions for the equation fl2.2ip . 

Eq. (12.211) is by no means the only equation that emerges from the set of relations 
(I2.17p . In fact, instead of singling out w = 5**^°'°^ we can chose other elements among 
the S''-*'-'-*, or (linear) combinations of them, and then systematically investigate what 
equations these choices satisfy by exploring the system of recurrence relations (I2.17p . For 
example, from (12.180 taking i = 0, j = — 1 and introducing the variable v = 1 — it 
is a simple exercise to obtain the following relation: 

^ pv-qv 

p — q-\-w — w = . [2.26) 

V 

Alternatively, adding the ^-shift of (l2T7al) to (]2.17dp we get 

pS^'''^ + qS^''^^ - /''^'^ + ^('•^■+^) = {p + q)S^''^^ + (5(''°) - 5 ''°^)§(°-^') , (2.24) 
and when taking in (12.240 z = 0, j = — 1 an easy calculation yields: 

^ pv + qv 

p + q + w — w = . (2.25) 

V 

Clearly, in (I2.25p . interchanging p and q and the ~-shift and the ^-shift should not make 
a difference, since the left-hand side is invariant under this change. Thus, the right-hand 
side must be invariant as well, leading to the relation: 

p (vv — vv] = q (vv — vv) . (2.26) 



Eq. (I2.26P is an integrable PAE in its own right for the quantity v, for which, by 
construction, we have an infinite family of solutions, namely given by 

v = l- 5^°'"^) = 1 - *c K-^ (1 + M)~^ r . (2.27) 

The PAE (12.260 for the variable v is identified as the lattice potential MKdV equation, 
which also occurred in [221 EZl, and is actually closely related to the lattice sine-Gordon 
equation of [18]. The relations (I2.23P and (I2.25P constitute a Miura transform between 
the lattice potential MKdV (12.260 and the he lattice potential KdV equation (I2.2ip . 

As another choice of dependent variables we can consider is the variable S^~^~^\ 
i.e. we can consider (12.170 for i = j = —1, leading to 



p 



+ = 1 - (l - (1 - ^(O'-i)) , 



and a similar relation for p replaced by q and the -shift replaced by the ^-shift. Using the 
fact that S^-^'O) = 5^° = 1-v and introducing the abbreviation z = S'^~^'~^'> ---- , 
the latter relations reduce to: 

p{z — i) = w , q{z — z)=vv . (2.28) 
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On the one hand, these two equations lead back to the equation fl2.26p by ehminating 
the variable z (considering in addition the ~- and ^-shifts of the two relations). On the 
other hand, by eliminating the variable v we obtain yet again a PAE, but now for z 
which reads 

^^-^^^-g = 4 , (2.29) 

which we identify with the Schwarzian lattice KdV equation, and is also referred to as 
cross-ratio equatior^. Here we have constructed iV-soliton solutions of eq. (I2.29p given 
explicitly by the formula 

z = 'cK^ (1 + M)-' K-'r-zo----, (2.30) 

P Q 

in which zq is an arbitrary constant. 

To summarise, we conclude that through the recurrence structure encoded in eqs. 
fl2.17p we obtain solutions of various different PAEs in one stroke. We will now proceed 
further by establishing a large parameter class of additional lattice equations, which also 
provide us information on the bilinear structure of the lattice systems. 

2.3. Bilinear Equations and the NQC Equation 

We will start by considering the r-function for the soliton solutions, given by (12.70 . using 
the relation ( 12.6ap . we can perform the following straightforward calculation: 

/ = det (i + m) =det{l+ [{pi + K)M + r'c] [pi + Ky^} 

= det{{pl + K) [1 + M + {pi + K)-^r*c\ {pi + Ky^} 
= det {(1 + M) [1 + (1 + M)-\pl + K)-^r'c] } 
= f det {1 + (1 + My\pl + Ky^r'c} 
from which, using also 



^ pl + K 

r = r 

pl-K 



we have 



y = 1 + *c (1 + M)~^ (pi - K)-^ r = 1 + (pi -KY^r , (2.31) 

where in the last step we have made use of the determinant relation (a special case of 
the famous Weinstein-Aronszajn formula): 

det (1 + X 7/^) = 1 + ■ a; , 

for arbitrary A^-component vectors a;, y (the suffix T denoting transposition). 

^ Eq. (|2.29p was first established as an integrable lattice equation in [5D], but was also studied in 
connection with discrete conformal function theory, [lOj . Interestingly, also this equation has appeared 
in the context of numerical analysis, in conection with the Pade tables in work by R. Cordellier, [13j . 
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The combination emerging on the right-hand side of eq. fl2.3ip is a new object which 
we need in the scheme, and it is natural to try and derive some equations for it along 
the lines of the derivations in subsection 2.1. Thus, more generally, let us introduce the 
function: 

V{a) = l-'c{a + Ky^u^""^ = 1 - (al + K)-^ r , (2.32) 
for any value of a parameter a G C. From f l2.3ip we immediately have 

j = ^ = V{-p) , (2.33) 

but we need further relations involving V{a) to derive closed-form equations for the 
r-function /. To do that we introduce some further objects, namely: 

u{a) = (l + M)~^(al + K)"V , (2.34a) 

'u{h) = 'c{hl + K)-\l + M)-\ (2.34b) 

5(a, b) = *c (61 + K)-\l + My\al + K)-^r . (2.34c) 

It can be shown that the latter object has the symmetry: 

S{a, h) = *c (61 + Ky^u{a) = 'u{b) {al + K)" V = S{b, a) , (2.35) 

in which a, 6 G C are arbitrary parameters. 

Following a similar derivation as the one leading to (12.111) and fl2.12p . derive the 
following relations: 

{pi - K)u{a) = i/(a)w(°) + {p- a)u{a) , (2.36a) 

{pi + K)u{a) = - V{a)u^°^ + {p + a)u{a) , (2.36b) 

Furthermore, by multiplying (12.360 from the left by the row vector *c {hi + K)^^ , show 
that fl236|) leads to 

l-{p + b) S{a, b) + {p- a) S{a, b) = V{a) V{b) . (2.37a) 

This and its companion equation 

l-{q + b) S{a, b) + {q- a) 5(a, b) = V{a) V{b) . (2.37b) 

form one of the basic relations for the rest of this paper. In fact, by setting a = —p, b = 
— g in fl2.37bl) we obtain 

1 + {P + q)S{-p, -q) = ^ = I , (2.38) 

which allows us to express the quantity S{p, q) in terms of the r-function. 

We will now now use both relations fl2.37ap and fl2.37bp . together with the the 
symmetry S{a, b) = S{b, a) , cf. (12.350 . to deduce now a partial difference equation for 
S{a, b) for any fixed a, b. In fact, from the identity: 

V{a) V{b) V{b) V{a) 



V{a) V{b) 



V{b) V{a) 
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by inserting fl2.37p and its counterpart, with p replaced by q and the ~-shift replaced by 
the ^-shift, as well as the relations with a and b interchanged, we obtain the following 
closed form equation for S{a, h): 

l-{p + h) g(a,&) + (p-a) g(a,&) ^ 1 - (g + a) S{a,h) + {q - h) S{a,h) 
1 - (g + 6) S{a, h) + {q- a) S{a, h) I - {p + a) S{a, b) + {p - b) S{a, b) 

In eq. (12.391) the parameters a and b are assumed fixed, and for each choice of them 
we have a quadrilateral PAE which is integrable in the sense of the multidimensional 
consistency property, explained in e.g. refs. [211 where p and q play the role of 
lattice parameters. We will elucidate the role of the parameters a, b in subsequent 
sections. It suffices here to observe that by fixing the special choice a = p, b = —p or 
a = q, b = —q respectively in (I2.37p . we obtain the relations 

V{p)V{-p) = l , V{q)V{-q) = l. (2.40) 

There are additional relations, connecting the object V{a) with the variable w defined 
in (I2.22p . which will play an important role in the proofs of the main results. Such 
relations can be obtained by eliminating the terms containing the product Ku{a) in 
the relations (I2.36al) and (]2.36bl) and their counterparts in the other lattice direction, 
and subsequently multiplying the the resulting combinations from the left by *c, using 
the fact that ^cK u{a) = 1 — V{a) , ^cuq = w. Thus, we arrive at the following list of 
relations: 

(2.41a) 
(2.41b) 

(2.41c) 
(2.41d) 

As a direct corollary the equality on the right-hand sides of eqs. fl2.41al) and (I2.41bl) . 
or equivalently (I2.41cl) and fl2.41dl) . actually provide us with another integrable lattice 
equation for the variable V{a). Once again, the parameter a plays a distinct role in this 
equation from the lattice parameters p and q, and in terms of the latter parameters the 
equation for V{a) is multidimensionally consistent. 

Remark: Note that closed-form equations for V{a) are obtained by equating the right- 
hand sides of f l2.41ap and (I2.41bl) . or equivalently the right-hand sides of fl2.41cp and 
fl2.41dp . Furthermore, breaking the covariance between the lattice directions by choosing 



p — q + id — w = {p + a) 
= {p-a) 

p + q + w — w = {p + a) 
= {p-a) 



V{a 
V{a 



V{a 
V{a 



V{a 
V{a 



V(a 



- {q + a 
-{q-a) 

+ {q-a 
+ {q + a 



V{a) 
V{a} 
V{a) 

V{a) 

V{a) 
V{a) 

V{a) 
V(a) 
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a = p in fl2.4ip we get for V{p) the following quadrilateral PAE 

Although this equation is not multidimensionally consistent in the strong sense 
(demanding consistency of the same equation in all lattice directions) it is 
multidimensionally consistent in a weaker sense (consistency between different equations 
on different sublattices). In fact, we can supplement (12.421) by a similar equation with 
a lattice variable h instead of m, and with lattice parameter q replaced by r. These two 
3-term relation are consistent-around-the-cube with a 4-term lattice equation for V{p) 
of the form arising from the right-hand sides of (]2.41ap . (I2.41bl) in the lattice directions 
associated with parameters q and r. 

Using now (12.331) to substitute the variable V{p), and similarly doing the same for 
V{q) in the analogous equation obtained by interchanging p and q and ~-shifts and ^- 
shifts, we obtain the following to bilinear partial difference equations for the r-function 
/, namely 

{p + q)ff+{p-q)f?= 2pff , (2.43a) 

{p + q)If + {q - P)fj = 2qff • (2.43b) 
Here the under-accents / and / denote lattice shifts in the opposite directions to / and 
/ respectively. It can be shown that these two 6-point equations are consistent on the 
two-dimensional lattice from an initial-value point of view, but we will not go into this 
here, cf. [23]. We just mention that by performing this computation explicitly, and 
eliminating intermediate values on vertices in the lattice one can derive from (12.431) the 
following 5-point lattice equation: 

{p - qfff- (P + qfll+ ^pqf = • (2.44) 
which is Hirota's discrete-time Toda equation, cf. [T8] . 

2.4. Lattice KdV and lattice MKdV 

We finish this section by presenting the actual lattice KdV, which is related to the lattice 
potential KdV equation (12.211) by considering differences of the variables w along the 
diagonals. These differences can be expressed in terms of the r-function, namely by 
setting a = p 01 a = —p in eqs. (12.411) . leading to the relations 

'E = p — q + w — w = {p — ^)^Jj ' (2.45a) 

T = p + q + w — w = {p + q)^-L , (2.45b) 

// 
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from which we have immediately 

ST=p2_g2 ^ H-S = T- f, (2.46) 

which leads to the lattice KdV equation, [18], in terms of either S or T, by eliminating 
one or the other variable using the first relation, i.e. 

E-E = {p'-q')(l-l) ^ T-T = (p^-g-^)|--^l . (2.47) 




What is not well known is that the resulting equation admits a scalar Lax pair of the 
form 

^=T(^+A(^, (2.48a) 
^=$ + E^, (2.48b) 

where A = A;^ — is the spectral parameter. In a similar way, from potential lattice 
MKdV equation fl2.26p . by considering ratios of the variable v over the diagonals in the 
lattice, we can obtain the following lattice equation for the variable W = v/v 

W ^ {pW -j) jp-qW) 
W (p-qW) (pW-q) ' 

The PAE (12.491) . which we identify with the (non-potential) lattice MKdV equation, 
arises as the compatibility condition of the following Lax pair 

^ =Wi^+ {q-pW)tp , (2.50a) 

p — qW p — qW 

where if) is a, scalar function and in which A is again a spectral parameter. For the 
variable W we have the following identifications in terms of the r-function: 

V = V{{)) = ^ ^ W = 3- , (2.51) 

/ fg 

in which g = T^f = det(l — iVf) is the shift in a direction with lattice parameter a = of 
the r-function. Using eqs. fl2.43p . which also hold for g, together with relations between 
/, g (which follow from (12.431) taking one or the other of the lattice parameters equal 
to zero) 

gf + fg = 2fg, (2.52a) 

gf + fg = 2fg. (2.52b) 

Finally, to obtain some explicit formulae for the r-function (12.71) we can use the 
properties of the Cauchy matrix M, in particular the fact that we can explicitly obtain 
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its determinant in factorised form. In fact, invoking the explicit determinantal formula 
for Cauchy determinants: 



det(A) 



(2.53) 



in which A is a Cauchy matrix with entries of the form: 

1 



hi Ij 



i,j = l,...,N 



where the ki, Ij, i, j = 1, . . . , N, are a collection of distinct parameters. Noting that the 
matrix M can be written as a Cauchy matrix of the form A with Ij = —kj multiplied 
from the left by a diagonal matrix with entries pi and from the right by a diagonal 
matrix with entries Cj, we have: 



det 



PiCj 



ki ~\~ kj 



PiCi I ■j-j / ki kj 
'^ki I , , \ ki -\- kj 



i<j 



(2.54) 



in which i,j can run over any subset of indices of {1, . . . , A^}. 

Furthermore, we have the following expansion formula for the determinant of a 
matrix of the form 1 + M 

N 

det (1 + M) = l + J2 l^^.^l + Yl 



i=l 



+ 5: 

i<j<k 



i<j 

Ma Mi 



Mu M, 



Mi 



Mj^i Mjj Mj^k 
Mk,i Mk,j Mk,k 



+ ■■■ + det(M) . 



(2.55) 



where as a consequence of fl2.54p we can compute all the terms in the expansion 02.551) 
explicitly, leading to a form similar to the celebrated iV-soliton form for the KdV 
equation given by Hirota, [T8] . 



3. From lattice KdV soliton solutions to Q3 solitons 

The main lattice equation arising from the scheme set up in section 2, i.e. equations 
that we coin "of KdV type" is the equation fl2.39p . which to our knowledge was first 
given in [22], cf. also [27], and which we can rewrite in affine linear form as follows: 



l + {p-a)S-{p + b)S 
l + {q-a)S-iq + b)S 



l + {p-b)S-{p + a)S 
l + {q-b)S-{q + a)S 



(3.1) 



Following [28] , we will refer to (13. ip as NQC equation. We begin this section by making 
precise the connection between the NQC equation and the equation Q35^q. This makes 
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some (mainly notational) ground- work for the subsequent statement and proof of the N- 
sohton solution for equation Q3, which is constructed on the basis of solutions of (13. ip . 
In [28] a full classification of all parameter subcases of (13. ip was given, In particular, it 
was already noted in the earlier papers that by limits on the parameters a, b (as either 
tend to zero or infinity) eq. (13. ip reduces to the other KdV type lattice equations, 
namely (I2.2ip . (I2.26P as well as (12.291) . Furthermore, as already remarked in [3], (13.10 
corresponds to the 6 = case of Q3. We will first make this connection more explicit, 
leading to a new parametrisation of Q3, and then show that in fact the A^-soliton 
solutions of (13. ip for different values of a, b together constitute a solution for the full 
case (i.e. 5 7^ 0) of Q3, which is the main statement in Theorem 1. 

3.1. Connection between QS^^q and the NQC equation 

We shall now relate eq. (13. ip to the special case of (ll.lcp with 6 = 0, which we indicate 
by Q3^^Q. We first introduce the following dependent variable 



,0 



u"„ ^ = Fn,m{a',b) {1 - {a + b)Sn,m{a,b)) , (3.2) 



in which 



where 



P= = (p^-a^)(p^-6'). ^ - a')(q' - b') . (3.4) 

This brings (13. ip in the form: 

where the lattice parameters have now become points p = (p, P) respectively q = (g, Q) 
on the (Jacobi) elliptic curve: 

p, q G r := {(x,X)|X2 = (x^ - a^){x^ - b^)}. (3.6) 

The full equation Q3 now reads 

P{uu + uu) — Q{uu + uu) = {p^ — q^) \ Cuu + uu) + ^pQ j (3-7) 

and this corresponds to the form of the Q3 equation (ll.lcp in the ABS list by the 
following relations between the original parameters p, q and the new parameters p, q: 



p'-b-" p_ {b'- a')p (&^ - a')q 

pz _ 1 — q'^ — a'^ 1 — 



while u = {iP' — a^)u , the latter being the dependent variable of the equation OLlcD . 
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3.2. N-Soliton structure for Q3 

The main result of the paper [7j was to give the A^-sohton solution for equation Q3, 
however without presenting there the full proof. Here we will present a constructive 
proof based on the machinery developed in section 2, which we think reveals some of 
the structures behind the Q3 equation and the solutions, and their connection to other 
lattice equations. 

Theorem 1: The N -soliton solution of Q3 (13.71) . which we denote by u^^'' = uH^m is 
given by the formula 

u'^^^ =AF{a,b)[l-{a + b)S{a,b)]+BF{a,-b)[l-{a-b)S{a,-b)] 

+ Cr(-a, b) [1 + (a - b)S{-a, b)] + Dr{-a, -b) [1 + (a + b)S{-a, -b)] 

Here S'(±a, ±6) = S'„,m(±a,, ±6) are the N -soliton solutions of the NQC equation (13.11) 
with parameters ±a, ±b as given in (12.351) . The function r{a, b) = Fn.m{0'^ b) is defined 
in (13.31) and A, B, C and D are constants subject to the single constraint: 

AD{a + bf - BC{a-bf = (3.10) 



Remark: The pseudo-linear structure of this solution, as an almost arbitrary linear 
combination of four different solution of the NQC equation with the parameters a, b 
changing signs, is remarkable, the various choices of signs apparently being connected 
to various choices of pairs of branch points of the elliptic curve (13. 4p of the lattice 
parameters. As we shall see in the unravelling of this solution, the a, b not only play 
the role of moduli of those curves, but as lattice parameters in their own right (w.r.t. 
additional "hidden" lattice directions). 

Proof: We will now go over the various steps needed to prove Theorem 1, all of which 
are based on the relations of the A^-soliton solutions to KdV type equations given in 
section 2. 

Step # 1. We first introduce a new associated dependent variable U^^^ = U^m given 
by: 

[/W =(a + b)Ar{a, b)V{a) V{b) + (a - b)Br{a, -b)V{a) V{-b) 

- (a - b)Cr{-a, b)V{-a) V{b) - (a + b)Dr{-a, -b)V{-a) V{-b), 

in which V{±a) = K,m(±a), V{±b) = K,m(±&) are defined by fl2:32|) and F{a,b) = 
Fn,m{.0'^b) is defined by (13.31) . We recall that the connection between the objects V{a), 
V{b) and 5'(a, b) is given by (I2.37al) . and that this relation can be covariantly extended 
to other lattice directions as for example in (12.37bp . Using these with the equations 
(I2.4ip we can derive the following set of important Miura type relations. 
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1 
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in (13.111) the following hold: 
p — q + w — w 



p + q + w — w 



U 
1 

u 



[Pu-Qu- (p^ 



1 
U 



Pu — Qu — {p^ — q^)u 

Pu — Qu — ijP' — (f)u 
Pu — Qu — {p^ — q^)u 



1 

U 



(3.12a) 
(3.12b) 

(3.12c) 
(3.12d) 



where w = w^m, defined in (I2.22p . is the N -soliton solution of the lattice potential KdV 
equation fl2.2ip . The relations f l3.12p hold for arbitrary coefficients A, B, C, D (i.e., 
without invoking the constraint (I3.10p ). 

Proof: This is by direct computation. In fact, to prove (I3.12ap . we consider 

Pu — Qu — {p^ — q^)u = 
= Ar{a, h) {p + a){p + b){l - {a + b)S{a, b)) + 

-{q + a)iq + 6)(1 - (a + 6)5(a, b) - (p'' - g2)(l - (a + 6)5(a, b)) 

= A{a + b)r{a, b) \{p -q)-{p + a){p + b)S{a, b) + {q + a){q + b)S{a, b) + (/ - q^)S{a, b) 



A{a + b)r{a,b) (p + a) 



(p + b)S{a,b)) + (p - a)Sia,b)) + 



+ ... 



-{q + a) (^1 - (g + b)S{a, b) + {q - a)S{a, b) 

= A{a + b)r{a, b) Up + a)V{a)V{b) - (g + a)V{a)V{b) \ + ... 
= A{a + b)r{a, b){p-q + w - w)V{a)V{b) + ... 
= {p — q + w — w)U , 

in which the dots in each hne on the right hand sides stand for similar terms with (a, b) 
replaced by (a, —6), (—a, 6), (—a, —6) and A replaced by B,C,D respectively. In the last 
steps we have made use of ( 12.37P and ( ]2.41ap respectively. The other relations (I3.12bp - 
( ]3.12dp are proven in a similar fashion. 

□ 



Step # 2. We now establish an important property of the objects defined in (13. 9p and 
(13. lip , which holds for arbitrary coefficients A, B, C, D (i.e., again without invoking 
the constraint (13.100 ). namely: 
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Lemma 2: The following identities hold for the N-soliton expressions: u = uh^Xi 



UU - P{u^ + n^) + {2p' - a' - h'')uu = ^ det(^) 



in which the 2x2 matrix A is given by 

A = 



(a + b)A {a - h)B 
-{a-b)C -{a + b)D 



(3.13a) 
(3.13b) 

(3.14) 



Proof: Again this is proven by direct computation. However, in order to make the 
emergence of the determinants more transparent we consider the following 2x2 matrices: 



M 



Pu - {p"^ - b'^)u , {p-b)U 

{p + b)U , -Pu + {p^ -b'^)u 

Qu-{q^-b^)u , {q-b)U 

{q + b)U , -Qu + {q^ -b'')u 

Evaluating the entries in these matrices we obtain 
Pu - (p2 - b^)u = 
P^ 



(3.15a) 
(3.15b) 



Ar{a,b) 



(a + b)S{a, b)) - {p' - b'){l - (a + b)S{a, b)) 



{p-a){p-b) 

Aria, b){p + b) [{p + a)(l - (a + b)S{a, b)) -{p-b){l-{a + b)S{a, b)) 
Aria, b){p + b)ia + b) [l - (p + a)Sia, b)) + (p - 6)5(a, 6)1 + . . . 
A{a + b)r{a, b){p + b)V{b)V{a) + ■■■ = ^p^-h^ r'^{a)Ar{b) , 



+ 



(3.16) 



in which again the dots in each hne on the right hand sides stand for similar terms with 
(a, fe) replaced by (a, —b), {—a,b), (—a, —b) and A replaced by B,C,D respectively. On 
the right-hand side of this expression we have introduced the vectors 



Via 



p^/^{a)V{a) 
p^/^i-a)Vi-a) 



r^{b) = {p'/\b)V{b),p'/\-b)V{-b)) 



in which the plane-wave factors p{a) are given by 

\p — a J \q — a ^ 
A similar computation as above yields 

-Pu + (p2 - b'^)u = ^p^ - 62 r'^{a)Ar{ 



(3.17) 



(3.18) 
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whilst U and U can be written as 

U = r'^{a)Ar{b) , U = r'^{a)Ar{b) . (3.19) 

Thus, we find that the matrix L in (13.151) can be written as 

^p2 _ ^2 r'^{^a)Ar{h) , { p - hy '^ {a)Ar{h) 
{p + b)r^{a)Ar{b) , ^p^ - b^r^ {a)Ar{b) 

and similarly for the matrix M. Using now the general determinantal identity 

det ^jVj^ = det {{yf ■ Xj)ij=i...^r) 

for any collection of r pairs of r-component column vectors Xi,y^ (the superindex T 
denoting transposition), we obtain the following result: 

det(lv) = (p2 - b^) det (A7{b) r^(a) + Ar{b) ?^(a)) 
= {p^ - b^) det{A) det (r{b) r^{a) + r{b) r^{c 

= det{(?W,.W)(^;[«j 

= - det (r(a),r(a)) det (r(6), r(6)) . 

It remains to compute the determinant of the matrix {r{a),r(a)) whose columns are 
the 2-component vectors r(a) and r{a). this is done as follows; 

det (r(a),?(a)) = p'/\a)p'/\-a)V{a) V{-a) - p'/\a)p^/^{-a)V{a) V{-a) 



p — a 
p + a 
p — a 
p + a . 



V{a)V{-a) 



p + a 



V{a)V{-a) 



p — a 

{p + a)S{a, —a) + (p + a)S{a, —a) 



+ 



p + a 
p — a 



1 — {p — a)S{—a, a) + {p — a)S{—a, a) 



2a 



where we have used ( ]2.37p for the choices b = —a and the signs of a,b reversed, as well as 
the fact that S{a, b) = S{b, a). Thus, putting everything together we obtain the result: 

det(X) = (/ - b') det{A) r^^^r.^, ■ (3.20) 

On the other hand a direct computation of the determinant gives: 

det(X) = -\Pu- (/ - b^)u\{Pu - (p2 _ 62)^ _ (p2 _ ^2-^^^ 

= -{P^- b^) \P{u^ + u^) - (2p' -a^- b^)uu - Uu] . (3.21) 

Comparing the two expressions for det(2v) from (13.201) and (I3.2ip we obtain the first 
equation in Lemma 2. 

□ 
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Step # 3. The last step is by combining the relations fl3.12p and fl3.13p as well as the 
potential lattice KdV equation (12.211) to assert that u = ul-j^m solves the Q3 equation. 
In fact, multiplying for instance ( I3.12bl) by ( ]3.12d[) and using ( I3.13b[) . where we identify 

52 



det(yi) 



16ab 



according to (13.101) . we obtain from the lattice potential KdV: 

p^ — q^ = {p + q + w — w){p ~ q + w — w) 

= Pu — Qu — [p^ — q^)u Pu — Qu — {p^ — q^)u 

{p^ - q 



Q{u^ + u)- (2g2 



b )uu 



4Q 



P^{uu) + Q'^{uu) + (p^ - q^fuu - PQiuu + uu) + 
- {P^ - q^)P{uu + uu) + {p^ - q^)Q{v^ + u ) 
P'^iuu + uu) + (Q^ - P^){uu) + (p^ - - PQiuu + uu) + 



{p^ - q^)P{uu + uu) + [p^ - q^)Q{u^ 



u 



-{2q' 



b^)uu + 



4Q 



P 



P{uu + uu) — Q{uu + uu) — (p^ — q^){uu + uu) 



+ 



+ {p" - 



2 2\~~ 

q )uu 



2 7 2 2 2 \ / ' 

2 +0 — p — q )uu + (p' 

where we have used the fact that p, q are on the elliptic curve (13.41) and the identification 

of UU via the result (13.131) of Lemma 2. From the last step, after some cancellation of 
terms, we obtain Q3, in the form (13.71) . for the function u 
proof of the theorem. 



u^nX^ which completes the 



Thus, we have established the explicit form of the A^-soliton solution ( 13. 9p of the full 
Q3 equation in the form (13. 7p in a constructive way, using basically all the ingredients 
of the lattice KdV type soliton structures of section 2. The parametrisation using the 
KdV lattice parameters g, in terms of which the underlying plane-wave factors (12.21) 
take a simple form, now involving some new fixed parameters a, fe, in terms of which we 
have an elliptic curve (13. 4p appearing in the equation, of which ±a, ±6 are the branch 
points. It turns out that it is natural now to covariantly extend the lattice to include 
lattice directions for which these parameters a, b are the lattice parameters, and thus 
to introduce extensions of the plane- wave factors of the form, [7j: 



Pn,m,a,p{ki 



p + hy fq + k 



p-kij \q~k 



a + ki 
a — kj 



b + ki 
b — ki 



p1 



(3.22) 
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This allows us then to use relations of the type fl2.33p and (12.381) to express the objects 
V{a) and S{a, b) in terms of a r-function / = fn,m,a,i3 as follows 

V{a) = ^ = ^ , l-{a + b)S{a,b)= = " ^ , (3-23) 

in which T^, Tb denote elementary shift in the four- dimensional lattice in directions 
associated with the parameters a, b (i.e. shifts by one unit in the corresponding lattice 
variables a, (3 respectively). Thus, as a consequence of Theorem 1, we get the following 

Corollary: The solutions (13. 9p takes on the following form in terms of the r-function 
f = fn.m,a,i3, whcn thc planc-wavc factors p are covariantly extended according to the 
formula (I3.22p . namely 

(N) Ai-r i,\ /n,m,Q-l,/9— 1 , r, ^ / i\ /n,m,a-l,/3+l 

'=Ar{a,b) — Y \-Br{a,-b) — V — 

Jn,m,a,l3 Jn,m,a,l3 
I /^i-r i,\ /n,m,Q:+l,/3— 1 . ^ / i\ fn,m,a+l,f3+l fn, n a\ 

+ Cr{-a,b) — '— \-Dr{-a,-b) — y — (3.24) 

Jn,m,a,l3 Jn,m,a,l3 

cf. whilst its adjoint function takes the following form: 
= A{a + b)r{a, b) ^"'"'"-y^"'"^-"'^'^ + B{a - b)r{a, -b) f^'^'-^ff-'^^'»''^+^ 

J n,m,a,l3 J n,m,a,(3 



r~il 7 N I- / iL\ fn,m,a+l,f3fn,m,a,B~l ,-,/ , i \ ^/ i\ fn,m,a+l,fifn,m,a,B+l 

C(a — b)F[—a, b) — ' — D[a + b)F[—a, —b) — ' — — 



J n,m,a,l3 J n,m,a,l3 

(3.25) 

As was indicated in [7], the verification of the solution of Q3 in terms of / can 
also be done by establishing a set of interlinked Hirota-Miwa equations in the four- 
dimensional space of the variables n, m, a, (3. An alternative approach to the iV-soliton 
solutions for ABS type equations is developed in a subsequent paper, [17], in which a 
formalism is developed in terms of Casorati determinant expressions for the r-functions. 

3.3. Associated biquadratic polynomials and Lax pair 

In this subsection we make some additional observations and connections related to the 
A^-soliton solution for Q3 presented above in Theorem 1. To this end it will be useful 
to introduce the following 



Qp,q(M, M, M, u) := P{uu + uu) - Q{uu + uu) - {p^ - q^) i^uu + uu) + -^pQ ^ 



(3.26) 



this quadrilateral expression is a polynomial of degree one in four variables with 
coefficients which depend upon the lattice parameters p, q G F (13. 6p . By introducing 
this polynomial the equation ()3.7p may be written conveniently as Qp,q(u, u, u, u) = 0. 
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Following the approach of ABS [3] we may associate to this polynomial a biquadratic 
expression which we denote "K 

?{p(M, u) := P{u^ + M^) - (2p2 - - b'^)uu + —. (3.27) 
This biquadratic is related to (13.261) in two ways, first, by the equations 

iq' - p'):K,{u,u) = Qp,q(Qp,q)- - (Qp,,)n(Qp,q)~ 

(3.28) 

where the subscripts u etc. denote partial derivatives. And, second, by the equation 

[Kp(n, u)^p{u, u) — ^q(u, u)'K^{u, u) = Qp,q(^^, u, u, u)Q* (j(u, u, u, u) (3.29) 

where Q* ^ denotes an associated quadrilateral given by: 

^ ^ ~ / ^ 5^ \ 

Q* q(M, u, u, u) := P{uu + uu) + Q{uu + nn) — (p^ + — — 6^) ( (mm + nn) — j ■ 

(3.30) 

Furthermore, the underlying curve associated to (13.261) . is determined by the polynomial 
in a single variable, 

r{u) = - 2%,{'K,U = - V - 6' . (3.31) 

The proof in subsection 3.2 relies on relations for a number of quantities that have 
appeared in passing, but which we will now explain somewhat further. In fact, the 
identities (I3.13P constitute factorisation properties for the aforementioned biquadratics, 
and can be written afl: 

'K^{u,u) = UU , 'K^{u,u) = UU . (3.32) 

As a consequence of the factorisation (13.321) and the identity (I3.29P it becomes clear 
that either the lattice equation Qp^q(n, n, n, u) = is satisfied, or the equation 
Q* q(u, u, M, m) = is. 

The matrices L and M appearing in the proof of Lemma 2, cf. fl3.15p . and 
whose determinants coincide with the biquadratics given above, are related to the Lax 
representation for Q3, which can be written as: 

:h T ( ~\A. ( Pu-ip'-k')u , -K \ 

= X,(m,m)0 = 7p(^^^^^^(^,_P^_^ -Pm + (^^^ - P)2 j ' ^'-^'"^ 

^ ( Qu-{q^-k^)u , -K \ / ,N 

= M,(M,M)0 = 7,(^^^^^^^,_p^_^ -gM+(g^-P)M j'^' (^-^^^^ 

^ Incidentally, by performing a natural continuum limit on one of the lattice variables, one can identify 
the object U with a derivative of u with respect to the emerging continuum variable. 
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in which we can set 7p = !Kp{u,u)~^/'^ , 7^ = 'Kq{u,u)~^^'^ . In [TTldn] it was explained 
how to obtain in general Lax representations for quadrilateral lattice equations which 
satisfy the multidimensional consistency property. In the case of the soliton solutions, 
choosing u as in fl3.9p . and using the corresponding factorisation of the biquadratics as in 
fl3.32p . with U given by (13. lip , one can explicitly identify the eigenvector cj) in terms of 
the ingredients of Theorem 1 and the objects of section 2. In fact, the construction of the 
Lax representation implies that in terms of the components of the vector cf> = {(pi, 02)"^ 
we have that the ratio 02/0i = T^u, where it is understood that all the variables are 
covariantly extended to a lattice involving a lattice direction associated with the spectral 
parameter k as lattice parameter. We will not pursue this identification in the present 
paper. 



4. Connection to the Backlund transformation 

In this section we study the recursive structure of the A^-soliton solutions for Q3 
presented in Theorem 1, under the application of the relevant Backlund transformation 
(BT). Thus, in Theorem 2 we give the precise correspondence between the A^- and 
soliton solutions, which, as we will see, will involve some redefinitions of the constants 
as we increase the soliton number. In particular, this is once again a manifestation of 
the multidimensional consistency of the lattice equation, by which we can interpret in a 
precise sense the lattice equation as its own BT. Furthermore, we highlight at the end 
of this section how the initial step of the Backlund chain, i.e. the seed solution defined 
as a fixed point of the BT, is embedded in the general formula of Theorem 1. 

4.I. Backlund transform from N- to N + 1-soliton solution 

The following theorem relates the A^-soliton solution given previously (13. 9p to the 
solution one would find by the iterative application of the Backlund transformation. 

Theorem 2: Let n*^^^ be as defined in l{3.9i) and letu^^^^^ he equal to u^^^^^ as defined in 
lis. 9\) . depending on additional parameters k^^i, Ctv+i and additional plane-wave factor 
Pn+i, cind where all hut the latter plane-wave factors pi, i = 1, . . . , N , as well as the 
discrete exponentials F(±a, ±6) replaced hy 

_ kN+i^ki Km+1 I ,x 

Pi = -r T:Pi , («^iV + l) , r±a, ±5= 77 —^77 —r:F[±a,±h), 

k^+i - ki [kN+i =F a) (/ctv+i =F 0) 

(4.1) 

respectively. Then u^^^ is related to u^'^'^^^ hy the Backlund transformation with 
Backlund parameter i^+i ^ that is the following equations hold 

Qp,....(«^^U^^^^^^^^^^^''"''^) = 0, Q,,,^,iu^^\u^^\u^^^'\%^''^'') = 0, (4.2) 
in which Qp^i,^_^_^ , Qq,t^^i are the quadrilateral expressions of the form given in (I3.26P 
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Remark: The shift pi — ^ corresponds to the covariant extension of the lattice 
variables in a direction given by the new parameter /cat+i, which equivalently can be 
described by the introduction of a new discrete lattice variable h associated with kN+i 
as a lattice parameter. This dependence can also be incorporated in the parameters 
Cj, as a hidden dependence on this variable. Remarkably, the coefficients A, B, C, D, 
which in principle could alter in the transition from the N- to + 1-soliton solution, 
in fact, remain unaltered in the incrementation of the soliton number. 
Proof: The proof of Theorem 2 can be broken down into a number of steps, which are 
all constructive, and rely once again on the machinery developed in section 2. In the 
first step we break down the + 1-soliton expression into components associated with 
the iV-soliton solution. In the second step we apply the BT to the A^-soliton solution 
and in the final step we compare the expressions obtained in the two previous steps. 



Step # 1. We first establish a recursive structure between the basic objects like 
Sn,m{ci,b) and Ki,m(o, 6) between the A^- and + 1-soliton solutions. This uses the 
breakdown of the Cauchy matrix as it occurs as the kernel (1 -|- A4")~^ in the various 
objects. 



Lemma: The following identity holds for the inverse of a (A^+1) x (A^+1) block-matrix 

-1 



A 


b 




d 



A~\l + \bc^A-^) 


-\A~^b 


-^c^A-^ 

s 





in which A is an invertible N x N matrix b and are A^-component vector column- 
and row-vector respectively , and d is a. nonzero scalar, where the scalar quantity s, 
given by 

s = d — c^A^^b , 

is assumed to be nonzero. 

Proof: The formula can be verified by direct multiplication. The matrix is invertible if 

s is nonzero. 

□ 

Let M^^^^^ be the (A^ + 1) x (A^ + 1) Cauchy matrix with parameters ki, . . . , /cat+i 
as defined in a similar way as in ( 12. ip . Applying the Lemma to compute the inverse of 
the matrix 1 + M^'^'^^^ , which can be decomposed as above by setting 

A = 1+mW , 6 = c;v+i(Wil+^)"V , c = pA.+ic^(Wl+^)-' , rf=l+^^^±i^, 
we have 



(1 + M) 1 + S ^CN+lPN+lU{kN+i) ^u{kN+i) 


-S ^CN+lU{kN+l) 


-S-'^PN+l^u{kN+l) 
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where u{-) and are given in fl2.34ap and (]2.34bp . and where s now takes the form 



^ = 1 + ^^^^ (1 - 2W,5(^)(Wi, Wi)) = , (4.3) 

where = det(l + M^^^^^) is the r-function of the (A^ + l)-sohton solution. 

Using now the definitions (I2.35P and (12.321) . we obtain the following recursion relations 
for the to + 1 soliton objects: 

(a + /C7v+i)(6+ fcTv+i) 
X {l-{b + kN+i)S'^''\kN+i,b)) , 

(4.4a) 

^(A^+i)(«) = v^(A^)(«) _ ^-'cN^iPN^i k^^^^s^N)^^^ y(^)(A:;v+i) , 

a + kN+i 

(4.4b) 

which hold for all a, h (not coinciding with — /cj, i = 1,...,A^). Setting h = /cat+i 
(assumed also not to coincide with any of the —ki, i = 1, . . . , N) we obtain as a corollary 
of these relations, the following identifications: 

^ l-{a + kN+i)S^^\a, k^+i) ^ V^^'KkN+i) , . 

' l-(a + A;^+i)^(^+i)(a,M V^(^+i)(W) ' 

which in particular implies that the ratio of 1 — (a + A;7v+i)5'(a, A^at+i) between its A^- 
and A^ + 1-soliton value, is independent of the parameter a. 

Step # 2. In this step we will apply the BT with Backlund parameter I = (/, L) G F 
to the A^-soliton solution defined in (13.91) . i.e. we want to solve the system of discrete 
Riccati equations for a new variable v 

Qp,i(w('^\w^^\t;,^J) =0 , Qq,((M(^),M(^),t;,t;) = 0, (4.6) 

cf. [6]. To solve this system we can reduce the problem by identifying two particular 
solutions, which are obtained by covariant extension of the known solution u^^\ By 
covariant extension we mean that all the plane wave factors p^, (z = 1, . . . , A^), as well 
as the discrete exponential factors Fri,m(=tfl; =t&); which enter in the A^-soliton formula 
(e.g. through the Cauchy matrix M = M^^\ and the vectors r), have a dependence 
on an additional lattice variable h associated with the lattice parameter /, such that 



and 



(±a, ±6) = r„,m(±a, ±b) 



{l-a){l-b) 



ri±a,±b) = {- ^- — -]r{±a,±b). (4.7b) 
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Here the elementary discrete shift in h is indicated by a ~, which gives meaning to the 
expressions u = Tiu, u = Tf^u, i.e. for the thus covariantly extended solution ^ we 
have 

—{N) [N) (N) (TV) 

7/ ^ 7/ 7/ ^ 7/ 

"'n,m,h "'n,m,h+l ' —n,m,h "'n,m,h-l ' 

As a consequence of the construction of Theorem 1, the following equations are satisfied: 



QM(n(^\n(^),n(^),i^'^^) = , Q,,i(n(^U(^),n(^),i^'^^) = , (4.8a) 
Qp,i(w(^\w^^\w(^\w^^^) = , Q,,t(w(^),S(^),M(^),i(^)) = , (4.8b) 

where ( ]4.8b|) holds because of the symmetry of the quadrilaterals. Having established 
the two solutions u^^'^ and u^^'^ of fl4.6p we can now find the general solution of that 
system in the interpolating form: 

V = '- , 4.9 

1 + r] 

in which rj is some function to be determined from the following coupled system of 
homogeneous linear equations: 

V _ Qp,((MW,MW,MW,i^'^^) 



7] _ Q,,i(M(^),MW,MW,i^'^^) 



(4.10a) 
(4.10b) 



V Q,,[(MW,MW,M(^),i(^)) ■ 

The compatibility of this system is equivalent to the compatibility of the BT 04.61) as 
a coupled system of discrete Riccati equations, which in turn is a consequence of the 
multidimensional consistency of the lattice equation. Using the explicit expressions for 
the quadrilaterals f l3.26p . eqs. (14.101) reduce to 

^ Pm(^) - (p2 _ /2)^(7V) _ ^^(JV) p + ; + ^ _ ^ 



V Pu(^'> - (p2 - P)u(^) - Lu^^^ p-l + w-w 
ri Qu(^) - (q^ - p)u(^) - Lu^N) q + l + w-w 



(4.11a) 
(4.11b) 



T] QuW - (g2 _ p)uW _ Lu^N) ^ 

where in the last step we have made use of eqs. fl3.12ap and the backward shifted version 
of (I3.12cp in the relevant lattice directions. Finally, using eqs. ( I2.41bl) and (I2.41cp with 
q —>■ I and w ^ w, etc. setting a = I, which lead to: 

p - I + w - w = [p - l)^ — , P + I + w- w = [p + I)—— , 

V{1) ^^^) 

we get 

7j fp + l\ V{l)V{l) fj fq + l\ V{l)V{l) 



V \P-lj V{l)V{l) V V{l)V{l) 
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which can be simultaneously integrated, yielding the following expression for the 
function 77: 

V = Vn,m = VO K„™(/)F„,^(/) . (4.12) 



^P — I J \(1 — I , 

Substituting the expression fl4.12p into (14.91) we have obtained the general solution of 
the system (14. 6p . It will be convenient to recast the solution (14. 9 p in the following form 



^W) = ^{N) _ ( 1(2/ + ^ _ y;)t/W , (4.13) 
\l+r]J L 



where we have inserted the expression for the difference which can be 

directly obtained from the backward-shifted relation of (I3.12cp by taking p = q = [ and 
identifying the shifts w and w (and similarly the shifts on u and U) with the shift w 
associated with the parameter [. 

Step # 3. We now set the Backlund parameter [ = ^at+i = {k^-^i, K^^i), which implies 
that the "-shift from now on is the lattice shift associated with the lattice parameter 
^AT+i. This enables us to reexpress the factors in (I4.4ap in the following way 

1 - (A:^+i + a)^(^)(a, fc^v+i) = V^'^H^tv+OZ^^H") , (4-14) 

which is an identity that holds for all a, and which follows from 

l + il~b)S}''\a,b)-{l + a)S^''\a,b) =V}''\a)V^''\b) , (4.15) 

by setting I = b = k^^i, which in turn is a covariantly extended version of (12.370 in 
the lattice direction associated with parameter [. Using the identity (14.140 in ( ]4.4ap and 
inserting the result into the expression (13. 9p for u^^~^^^ we find 

^{N+D ^ ^ir^a, b) [1 - (a + 6)5(^+^)(a, b)] + ■ ■ ■ 



AF{a,b) 



(1 - {a + b)S^^\a^b)) - + 

[a + kN+i)[b+ kN+i) 



= - -^{V^''\k^+i)fpN+iCN+iS-' Aia + b)ria, b)V}''\a)V}''\b) + ■■■ 
Kn+i l 

where again the dots stand for the remaining terms with coefficients B, C, D and with 

(a, b) replaced by (a, —6), (—a, 6), (—a, —6) respectively, leading to the following formula 

for the + 1-soliton solution: 

^(A^+i) = _ -}—(V^^)(k^^,)Ypj,^,cj,^,s~'U^^^ . (4.16) 

On the other hand, setting [ = i-N+i, and choosing the constant 770 = CN+i/{2kN+i) in 
(I4.12p . we obtain the identification 

l + r] = s, (4.17) 
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with s given in fl4.3p . using also (14.141) in the case of a = /cat+i. Inserting (14.171) into the 
solution V of the BT, (14.131) . with [ = ^at+i, and using the relation 

2/CAr+i + W — W = ZKn+I- 



which follows from the backward shift of ( 12.41cP setting p = q = a = k^+i, we can 
make the identification v = m^^^^-* between (I4.13p and the forward shift of (14.161) . This 
establishes the relation (14.21) . and hence completes the proof of theorem 2. 
■ 

Theorem 2 establishes the precise connection between the structure of the A^-soliton 
solution as given by the Cauchy matrix approach, and the way to generate a soliton 
hierarchy through Backlund transforms. What we conclude is that these two approaches 
coincide up to a subtle identification of the relevant constants in the solution. Since, as 
was remarked in the Corollary of section 3 that the soliton solutions of Q3 really live in 
an extended four-dimensional lattice, the precise identification of those constants is of 
interest, since they contain possibly the additional lattice directions. In fact, in [6] we 
established the first soliton type solutions for Q4 through the Backlund approach, and 
it is of interest to see how that approach connects to a (yet unknown) representation of 
multi-soliton solutions in terms of a scheme similar to the one set up in this paper for 
Q3. 

4-2. Connection to the Backlund transformation fixed-point 

In [6] we introduced a method for the construction of an elementary solution (of 
a multidimensionally consistent lattice equation) which is suitable as a seed for the 
subsequent construction of soliton solutions by iterative application of the Backlund 
transformation. This method is based on an idea which can actually be traced back 
to Weiss in [30], cf [8] for a further discussion of this point. We will now describe the 
connection between the solution which arises when that method is applied to construct 
a seed solution for Q3 (13. 7p and the 0-soliton solution found by substitution of = 
into the principal object studied here (13. 9p . Note that when we substitute iV = in (13.90 
the functions S{±a, ±6) = ,„(±a, ±6) appearing in the solution all vanish identically. 

Following [6] we consider a solution of (13. 7p ue = {ue)n,rn which is related to itself 
by the Backlund transformation of (13. 7p with Backlund parameter t G F (recall F was 
defined in (13.60 ). I.e., uq satisfies the system 

%,t{ue-, Ue, Ue, ue) = 0, Qq,tiue, ue, ue, ue) = 0. (4-18) 

These are coupled biquadratic equations for ue which are parameter deformations of the 
biquadratic introduced earlier (I3.27p . In this case we need to introduce a limit on the 
curve, F 3 t = (t,T) — > oo+ in which t — > oo and T ^ t^ - - + 0{t-^), 
then the identity 

lim Qp^i{u,u,u,u) = —0-Cp{u,u), (4.19) 
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reveals the sense in which the one is a deformation of the other. Recall [B] that with this 
special choice of t the system (14.181) defines a non-germinating seed solution (or singular 
solution in [1]). 

To solve simultaneously the biquadratic equations (14.181) (for general t) we introduce 
parameters pe and qe defined in terms of p and q by the quadratic equations 

Pe + 1/Pe = 2 qe + l/qe = 2 (4.20) 

The (canonical) solution of (14.181) may then be written 

ue = {ue)n,m = Aep'^eq^ + Bep^'^qe"', (4.21) 
where the coefficients Ag and Bq are constants subject to the single constraint 

" 16r(2t2 _ 2T - a2 _ • (4.22) 

We are now in a position to describe the connection between the solution found 
as a fixed-point of the Backlund transformation (]4.2ip and the solution u^^\ i.e., the 
solution obtained by setting = in (13.91) . Whereas that solution contains four terms 
with coefficients A, B, C and -D, the solution (I4.2ip contains only two terms, so it is 
clear the solutions do not coincide. Actually what we find is that making the particular 
choice t = (0, ah) the solution (^^21]) becomes with A = Ag, D = Bq and 5 = C = 0, 
whilst making the choice t = (0, —ah) (I4.2ip becomes m'-'^^ with A = D = {)^ B = Ag and 
C = Bg. In both cases the constraint (I4.22p becomes (I3.10p . 

5. A^-soliton solutions for the degenerate sub cases of Q3 

The coalescence scheme illustrated in figure [3] will be used in this section to construct 
A^-soliton solutions for the equations Q2, Ql, H3, H2 and HI by degeneration from 
the iV-soliton solution we have given for Q3. We begin by detailing how the equations 
themselves are found by degeneration, then we show how to find the new solutions. 

Q3l -[021 -[qT 



H3 



H2 



HI 



Figure 3: Coalescence scheme employed to construct A^-soliton solutions for the 
degenerate sub cases of equation Q3. 
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5.1. The degenerations from Q3 

We degenerate from Q3 as parametrised in ( 13. 7p which differs from the parametrisation 
given originally by ABS [3] (which we reproduce in (11. Id) ). The degenerate equations 
consequently appear in a parametrisation different than in the lists (11.11) and (II. 2p . 
Importantly, throughout the scheme depicted in figure [3] there is no limit taken on 
the parameters p and q present in (13. 7p . so each equation emerges in terms of these 
common lattice parameters (which later will be seen to occur naturally in all of the 
7V-soliton solutions). These parameters are simply related to the parameters p and q of 
the equations listed in (11.11) and (II. 2p by the following associations 

p^ — a-^ F q'^ — a'' Q 

Q2,Q1: p=—^ ^, ^ = —2 J' 

p'^ — a-^ q-^ — 

H2, HI: p= -p\ q = -q^. (5. Id) 

When written in terms of these common parameters p and q, the equations Q3 and H3 
involve also upper-case parameters P and Q which are defined in terms of p and q by 
an algebraic relation. For Q3 this relation was introduced previously (13.41) and is given 
again in (15. lap , for H3 there is a different algebraic relation which is given in (IS.lcp . 

It is by limits on the parameters a and b and on the dependent variable u that 
the degenerations appearing in figure [3] are achieved. Specifically, the following list of 
substitutions result in the indicated degenerations in the limit e — > 0: 

> = a(l - 2e), M ^ A + 1 + (1 + 2u)e^ , (5.2a) 

6^ 1 

u ^ — + -u, (5.2b) 



Q3 - 


^Q2 


Q2 - 


-^Ql 


Q3 - 


H3 


Q2 - 


H2 


Ql - 


HI 


H3 - 


H2 


H2 - 


HI 



4e2 



e 



1 ^ . ^ 

6 = -, u-^ M, (5.2c) 
2 

a = -, u ^ — h e^M, (5. 2d) 
e 4 

a = -, M — > e5u, (5.2e) 

e 

1 ^ / e4 - 



c2 ■ 



u V-Se ( 1 + — M ) , (5.2f) 



u 



- + -u. (5.2g) 



e 

Note that for these degenerations we assume the parameter 5 appearing in the equations 
Q3, Ql and H3 is nonzero. 
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5.2. N -soliton solutions of Q2 and Ql 

We degenerate from the previously established A^-soliton solution for Q3 ( 13. 9p . To begin 
we will consider in detail the degeneration from this solution to the A^-soliton solution of 
Q2. We are led by the requirement that the parameter b which appears in this solution 
should, according to fl5.2ap . be replaced with b = a{l — 2e). Making this substitution 
and expanding the result in powers of e results in a lengthy expansion of (13. 9p which 
we break down into expansions of the simpler component parts. First we consider the 
function F{a,b) defined in (13.31) which has the following expansions depending on its 
choice of argument 

ir(a, b) p(a)(l + ei + e\e/2 + x)) + ^(e^), 
r{a,-b)^ l-e^ + e\e/'2-x) + 0ie'), 
r{-a,b)^ l + e^ + e\e/^ + x) + 0{e'), 
F{-a, -b) p(-a)(l -ei + e\e /'^ " x)) + 0(6^), 

in which we have introduced the new functions 



(5.3) 



y^a"^ — p'^Y (a^ — g^)"^ 



(5.4) 



and the function p(a) = pn,m{,o) coincides with the one defined previously (13.17p . In a 
similar fashion we find the expansion of terms involving S{a, b) defined in (12.341) to be 
the following 

I - {a + b)S{a,b) — > 1 - 2a^(a, a) + 0(e), 

I- {a-b)S{a,-b) — > I - e2aS{a, -a) + e^Aa^ Z{a, -a) + 0{e^), 

l + {a-b)S{-a,b) — > l + e2aS{-a,a) + e^Aa'^Z{-a,a) + 0{€^), 

1 + (a + b)S{-a, -b) — > 1 + 2a5(-a, -a) + 0(e), 

in which we have introduced the new quantity 

Z(a, b) = *c {bl + K)-^u{a) = *c {bl + K)~^ (1 + M)'^ (al + K)-^ r. (5.6) 

Now, although the expansions of these components of (13. 9p are fixed, by choosing the 
dependence of the constants A, B, C and D in (13.91) on the small parameter e we can 
exert some control over its overall expansion. In order that it be of the required form, 

namely u^^'' — > — - ( — h 1 + (1 + 2M*^^))e ) as listed in (I5.2ap . we make the following 
4a^ \e / 

choices for these constants 
A 

^^3^(- + l-^o + ((3 + eo')/2 

(5.7) 

O^ — (- + l+eo + ((3 + eo')/2 
5 



5 I 










-1 

8a2 ' 




5 , 
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Here the four constants A, B, C and D (constrained by fl3.10p ) which are present in 
the Q3 A^-sohton solution (13.91) have been replaced with three modified constants A, 
D and (with no constraint), in the degeneration to the Q2 A^-soliton solution. The 
number of these constants minus the number of constraints is therefore preserved in the 
degeneration, which is strong evidence that the solution found by degeneration is the 
most general one obtainable by this method. Finally then, the solution of Q2 which 
arises in the limit e — > as a consequence of the substitutions (15.31) . (15. 5p and (15. 7p 
into (13.91) is found to be 

=7((^ + ^o)' + 1) + a(e + Co)S{-a, a) + {Z{a, -a) + Z{-a, a)) + 
4 ^ ^ (5.8) 

AD + -Ap{a){l - 2aS{a, a)) + -Dp(-a)(l + 2aS{-a, -a)). 

Here ^o, A and D are the aforementioned constants which may be chosen arbitrarily, 
^ = ^n,m and p(a) = Pn,m{ci) are defined in (15.41) and (I3.17p . and S{a, h) = Sn,m{ci, h) and 
Z{a,b) = Zn,rn{ci,h) are defined in (12.340 and (15.61) respectively. 

We now consider the degeneration Q2 — > Ql which is somewhat simpler than the 
Q3 — > Q2 degeneration considered above. To achieve the required limit fl5.2bp of the 

Q2 A^-soliton solution (15. 8p . namely that m*^^) — > — - + -u^^\ we need only substitute 

Ae'^ e 

the constants appearing in the solution as follows 

D (5.9) 



e 



2B 
e 



The desired limit results from substitution of (15. 9p into (15.80 . provided the modified 
constants A, B, D and .^o are chosen to satisfy the single constraint 

AD + \b' = ^. (5.10) 

Notice that again the number of constants minus constraints is preserved in the Q2 — > 
Ql degeneration. The A^-soliton solution of Ql which emerges is 

«W = Ap{a){l-2aS{a, a)) + B{^+^o+2aS{-a, a))+Dp{-a){l+2aS{-a, -a)). (5.11) 

As a small side development here we now consider the further limit a — > in the 
Ql A^-soliton solution (15. lip . There is some subtlety in this limit and the solution of 
Ql which emerges inspires some useful additional observations. Performing the a — > 
limit naively would also change the equation (in fact send it to Ql5=o) because a appears 
in its parametrisation fl5.1bp . However, the substitution 



a = e, u —>■ 1 + e^u, 



(5.12) 
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preserves the full equation as e — > 0. This leads to a reparametrisation of Ql in which 
we identify the parameters p and q present in fll.lal) simply as 

P=4, ^ = 4, (5.13) 

which come to replace the associations in (iS.lbl) . Now the same substitution, a = e, 
yields the following small-e expansion for the component parts of the solution (15. lip . 

p{a)^ l + eiy + e% + Oie^), 



p{-a)^ 1 -ez/ + e2| + 0(e3), 

-eu + 0{e'), 
S{a, a) - 2e5(-i-2) + q^^2^^ 

S{-a, -a) + 2e5("i'-2) + 0{e^), 



(5.14) 



where we have introduced the new function 



2 2 

-n H — m. (5.15) 
P Q 



If we also make the choice 



^+j{A + uo), 

l-- + ^((l-5)z/o-2MM), 

5-1 \ (5-16) 
-— + -A(2z/i/i.o-l), 

^0 ^ 1 - ez/Q + —vl, 

for the constants appearing in (15. lip then the substitution of fl5.14p yields the desired 
limit of the solution u^^^ — > 1 + e'^u^^\ The resulting solution of Ql is as follows 

= 5 Q(i^ + uof + z/o)5(-i'-i) + 25(-i'-2)^ + A{u + v,- 2S^-'^-^'^), (5.17) 

where the constants vq, vi and A are arbitrary, v = Vn.m is defined in fl5.15p and 
^(^.i) = S'^;^) is defined in (ESED. 

The main point we wish to make about the Ql solution fl5.17p is that it generalizes 
the solution previously given f l2.30p for the lattice Schwarzian KdV equation, i.e., the 
equation Ql5=o. Specifically this can be seen as an extension of that solution to the case 
(5 7^ in that it reduces to (12.30P if we take 5 = and A = —1/2. 

5. 3. N -soliton solution of H3 

To find the A^-soliton solution for H3 we degenerate from the Q3 A^-soliton solution 
(13. 9p led now by the requirement that, according to (I5.2cp . we choose b = —. As before 
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we give small-e expansions for tlie component parts of the Q3 A^-sohton solution fl3.9l) 
which result from making this substitution for b. We find that 

r{a,b) ^? + 0(e2), 

F{-a, b) — > (-1)"+"^^-! + 0(e2), 

r(-a,-6) ^ ?9-i + 0(e2), 

l-{a + b)S{a,b) — > V{a) + 0{e^), 

l-{a-b)S{a,-b) — > V{a) + 0{e^), 

l + la-b)sl-a,b) — > r(-a) + 0(e2), 

l + {a + b)S{-a,-b) — ^ vl-a) + ole% 

where we have introduced the new function 

^ = ^n,n. = (^] (^] , (5.19) 



(5.18) 



a — p / \a — q 

which involves parameters P and Q which are related to p and q by fl5.1cp . and where 
V{a) = Vn,m{a) is defined in ( 12:321) . 

Substituting the expressions (15.181) into (13.91) whilst choosing the constants in that 
solution to be 

A^e^^A, B^e^y^B, 

we find u^^^ — > e^—^u^^'' as required for the Q3 — >H3 limit given in (I5.2cp . Thus we 
find the A^-soliton solution of H3 to be 



u 



(TV) 



{A + B)W (a) + ((-1)"+™C + D)'d~W{-a). (5.21) 



Here the constants A, B, C and D are subject to the single constraint 

AD-BC = — , 
4a 

which follows by substitution of 6 = — and (15.201) into (13.101) . The functions d = 'dn,m 
and V{a) = K,m(a) are defined in fICTD and ^Ml- 

5.4- N -soliton solutions of H2 and HI 

To find the H2 A^-soliton solution we choose to degenerate from the Q2 A^-soliton solution 
(15. 8p led by the requirement that we substitute a = - (cf ( ]5.2d[) ). (Observe from figure 
[3] that we could choose to degenerate from the H3 A^-soliton solution (15.211) to the H2 A^- 
soliton solution, the two paths actually lead to the same result.) Making the substitution 
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a = - into the component parts of (15.81) yields the following small-e expansions: 
e 

p{a) + eC + + 0{^)), (5.22) 

pi-a) — (-1)-+-(1 -eC + + 0{e')), 

in which we have introduced a new function (, 

C = Cn,m = 2np + 2mg, (5.23) 

and 



(5.24) 



aSia, a) eS^"'") - 2e2^(o,i) + ^(^3)^ 
aSi-a, -a) eS^"'") + 2e^S^^'''> + ©(e^), 
a2(Z(-a, a) + Z{a, -a)) 2e^S^'>^'^ + 0{e^). 

Substituting (I5.22p and (I5.24p into (15. 8p combined with the following choice for the 
constants 

A^ A(e + CieV2), 

A(-e + CieV2), (5-25) 
^0 eCo, 

results in an expansion for (15.81) of the required form 
with the new H2 A^-soliton solution which results being 



= ^(C + Co)' - (C + Co)^^°'°^ + 25(°'i) - ^2 + (-1)-+-A(C + Ci - 25^°'°^ (5.26) 



,W = 

4' 



where the constants A, (q and Ci are arbitrary (and unrelated), ( = (n.m is defined in 
fl5:23|l and S^''^^ = sLi'i! is defined in (KEc^ . 

We remark that the solution (I5.26P in the case A = can be transformed to the 
solution given previously for the equation Ql (I5.17p with A = and 5 = 1 by the simple 
(self inverse) transformation 

p,q,h...kN^l/p,l/q,l/h...l/kN, pi^(-l)"+>i. (5.27) 

The connection between these solutions reflects a kind of duality between the equations 
Ql5=i and H2 which was found previously in [9]. 

To find the A^-soliton solution for the equation HI we choose to degenerate from the 
A^-soliton solution of Ql (15.111) . According to (I5.2ep we should make the substitution 

a = - into (15.111) . conveniently we have already expanded the component parts of this 
e 

solution in powers of e because they appeared in our consideration of the Q2 — > H2 
degeneration detailed above in (15.221) and (I5.24p . The required limit of the solution, 
which according to (I5.2ep is u^^^ — > e5u^'^\ is achieved by choosing the constants 
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appearing in (15.111) as follows: 



B 



D 



A 




+ Cie), 



(5.28) 



The resulting A^-soliton solution of HI reads 



u 



W = i?(C + Co-2^(°'0)) + ( 



(5.29) 



where Co? Ci? ^ ^^id B, subject to the single constraint 




are otherwise arbitrary constants, ( = (n,m and S^'^'^^ = Si 



are defined in (15.231) and 



( ]2.8cp respectively. 

5.5. Equations A2 and Al 

In the above we have constructed solutions by degeneration following the coalescence 
diagram of figure [31 This diagram does not include the equations A2 and Al and we 
have thus far not given explicitly their A^-soliton solutions. However these equations 
are related to Q35=o and Ql respectively by straightforward gauge transformation, so 
the solutions we have given for those equations may be transformed to solutions for the 
equations A2 and Al. 

6. Concluding remarks 

In this paper we have reviewed the construction of A^-soliton solution for integrable 
quadrilateral lattice equations of KdV type dating back to the early 1980s, cf. [221 127] . 
as well as constructed A^-soliton solutions for the majority of equations in the ABS list. 
We have concentrated particular on the case of Q3, which by degeneration yields all the 
other ABS equations (except Q4) by limits on the parameters. The A^-soliton solution of 
Q3 is particularly interesting as it is most conveniently described in a four- dimensional 
lattice, associated with four lattice parameters, two of which are the lattice parameters 
of the equation supplemented by two further parameters acting as branch points of an 
elliptic curve. The emergence of this elliptic curve is rather mysterious at the level of 
Q3, which, unlike Q4, does not really warrant an elliptic parametrisation. Nonetheless, 
this curve naturally arises through what we consider to be the universal parametrisation 
of all ABS equations (except perhaps Q4, which remains to be investigated), and which 
allows us to treat all equations in the list on the same footing (unlike the original 
parametrisation from [3] where the parameters of the different equations in the list do 
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not seem to be directly linked to each other). Thus, the A^-soliton solution of Q3 can 
be written as a linear combination of four terms each of which contains as essential 
ingredient the A^-soliton solution of the so-called NQC equation of [22] with different 
values of the branch point parameters which enter in that equation. The A^-soliton 
solutions of the other equations in the list, namely Q2, Ql, H3, H2, HI, follow by 
degeneration, and were derived in explicit form from the solutions of Q3. The present 
paper concentrated on an approach, using a Cauchy matrix representation of the basic 
objects, and which essentially was developed, as direct linearization approach in the 
paper of the early 1980s. As a remarkable upshot a novel Miura transform between Q3 
(almost at the top of the ABS list) and HI (at the bottom of the list) played a crucial 
role in the mechanism behind the solutions. In the subsequent paper [17] an alternative 
representation of the soliton solutions, in terms of Casorati determinants, is given, based 
on the bilinear forms of the ABS equations. 

We have not touched in this paper on continuum limits of the equations and 
the hierarchies of continuous equations associated with the lattice systems. This can 
obviously be done without problem. The most direct way of introducing continuum 
analogues is by simply replacing the discrete plane- wave factors pi, which as functions of 
the discrete variables n,m was given in (12.21) by exponentials, making the replacements: 



or by simply including the exponentials in the pi together with the discrete exponential 
factors (invoking once again multidimensional consistency). We leave the derivation of 
the corresponding PDEs as an exercise to the reader. 
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